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We consider lepton flavor violations (LFV) mediated by photino as a result 
^h. of the nondiagonal slepton mass matrices in general supersymmetric models. 



Using the experimental upper bounds on I — > I +7 and fj, + Ti —> e + Ti 



^ . as constraints on the flavor changing slepton mass insertions, we predict the 

possible ranges of the upper limits on the branching ratios of other LFV 
processes such as I — > 31' , r — ► l^+n (01 n, p°, 0), Z° — ► ijj^i, and muonium 
— > antimuonium conversion. Most of these decays are expected to occur with 
small branching ratios far below the current or future experimental search 
limits. We also derive constraints on the flavor conserving mass insertions 
from the anomalous magnetic moments of leptons. 
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I. INTRODUCTION 



At present, the minimal supersymmetric standard model (MSSM) is widely considered as 
the leading candidate for the physics beyond the standard model (SM) |L]. It can solve the 
gauge hierarchy problem by supersymmetrizing the SM (with an addtional Higgs doublet). 
The boson loop contribution to Higgs mass is cancelled by the fermion loop contribution, 
the latter of which comes in with opposite sign to the former. In doing so, the particle 
spectrum of the theory becomes doubled compared to that of the SM. One expects a lot of 
new scalar particles (superpartners of the SM fermions) and new fermions (superpartners of 
the SM gauge bosons and Higgs). These new particles should have masses around 0(100) 
GeV —0(1) TeV in order to solve the gauge hierarchy problem in terms of softly broken 
super symmetry. Nice features of the supersymmetric theories are their calculability using 
perturbation theory, and the decoupling nature of the loop effects of new (super)particles 
on various electroweak observables, except for the dangerous supersymmetric (SUSY) flavor 
changing neutral current (FCNC) and SUSY CP problems (which will be dicussed shortly 
in more detail). Therefore the successful predictions of the SM do not change very much 
even if we have doubled spectrum of particles in SUSY theories modulo SUSY FCNC and 
CP problems. 

However, in generic supersymmetric (SUSY) models, one has to pay for this extra sym- 
metry. First of all, the lepton family numbers (Lj =e /J T ) and the baryon number (B) are no 
longer conserved as in SM. One can write down renormalizable superpotential which violates 
the Li and B numbers and leads to too fast proton decay in conflict with the observation. 
Secondly, the soft mass terms for sfermions can lead to large FCNC unless certain condi- 
tions are met. In most phenomenological SUSY models, one solves the first problem by 
assuming R— parity conservation by hand. The second problem (SUSY FCNC) is solved by 
assuming that either (i) the sfermion mass matrices are proportional to the unit matrix in 
the flavor space @ || , (ii) the sfermion mass matrices are proportional to the corresponding 
fermion mass matrices so that both can be diagonalized simultaneously 01], or (iii) assum- 
ing that the first two generation sfermions are highly degenerate and very heavy (gtrsim50 
TeV >> M S usy so that they basically decouple) ||. In the minimal SUGRA models with 
the flat Kahler metric at Mpi anc k scale, the first condition can be met, namely the squarks, 
sleptons and Higgs are all degenerate at the Planck scale. However, when one evolves the 
sfermion mass parameters to the electroweak scale using renormalization group (RG), the 
off-diagonal elements of the squark mass matrices are induced in a calculable manner, al- 
though there is no LFV induced at low energy. Moreover, the condition of the flat Kahler 
metric is a strong assumption which may not be true in general. For example, SUGRA ra- 
diative corrections to the boundary conditions at M Planck scale induce generically 0(~ 10%) 
off-diagonal sfermion mass matrix elements |J. Therefore, one can imagine certain amount 
of nondiagonal sfermion mass matrix elements at the electroweak scale in general. In the 
lepton sector of the minimal SUGRA model, there is no LFV as in the SM since neutrinos 
are massless. As an example, consider the minimal SUGRA model with the flat Kahler 
metric. Then the scalar masses are universal (being wig) at the Planck scale, whereas at the 
weak scale the sfermion masses change as 

(m^LL^fi = m z ) = m d m\ + ml + c d m u m\, (1) 
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as a result of renormalization ||. Because of the last term containing m u mj^, it is not 
possible to diagonalize m d m} d and {m 2 j) LL simultaneously. This leads to the flavor changing 
gluino-quark-squark vertices, which can contribute to various low energy FCNC processes. 
Also this {LL) mixing induces the LR and RR mixings in the minimal SUGRA models. For 
sleptons, on the other hand, we have 

(mf) LL ([i = mz) = mimj + m 2 + cim v m\, (2) 

and massless neutrinos (namely, absence of righthanded neutrinos) imply that there is no 
generation mixing in the slepton mass matrix. Since LR and RR transitions are proportinal 
to the {LL) mixing, there will be no lepton family number violation in the minimal SUGRA 
models. However, if there is LFV at high energy scale (e.g., in Supersymmetric Grand 
Unification Theories (SUSY GUT)) [0 or if right-handed neutrinos are included in the 
SUGRA ||, there can be generic LFV at electroweak scale. 

In view of this, it is important to see how large deviation from the above conditions 
(i) and (ii) are allowed in the general SUSY models by the various FCNC processes at 
low energy. Such studies have been done previously by several authors already, mainly on 
the gluino-mediated FCNC in the quark sector and the photino-mediated — > 1^ + 7 [§J • 
Basically deviations between the first and the second families should be very small. In terms 
of a dimensionless parameter defined as 

{6 l ) AB = A'Wmf, (3) 

where m~ is a suitable average of the slepton masses, the condition that deviations be- 
tween the first and the second families should be very small can be represented as following 
constraints 0: 

(S[ 2 )ll = 0(KT 3 ), and {5 l 12 ) LR = 0(1(T 6 ). (4) 

On the other hand, the deviations involving the third family are more loosely constrained : 

(4(23)Ul = O(1-10), and {6 1 13{23) ) LR = 0(1(T 2 ). (5) 

This is in part due to the less precise experimental informations on various FCNC processes 
involving the third family. But there are many interesting possibilities for which one can 
treat the third famly in a different manner from the first two families. In such theories, 
one may expect larger deviations from the degeneracy in general, and thus expect FCNC 
processes with branching ratios that may be accessible in the near future. 

In this work, we mainly concentrate on the photino-mediated FCNC processes in the 
lepton sector, which are almost parallel to the work by Gabbiani et al. [JJ|. Namely, we 
assume that the slepton mass matrices are not diagonal in the basis where U — Ij—'y vertex 
is flavor diagonal. In order to simplify the analysis, we make an assumption that the lightest 
superparticle (LSP) is a photino (7), and other neutralinos are fairly massive so that their 
effects are negligible compared to the LSP effects considered in this work. Finally, we 
assume that the off-diagonal mass matrix elements of sleptons are small enough that the 
mass insertion approximations are applicable. All of these assumptions are the same as 
Ref. 0, except for the photino mediated LFV instead of gluino mediated FCNC. In the 
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case of glino-mediated FCNC, the neutralino effects will be generically suppressed by c^/as, 
so that one can safely ignore the neutrlino-mediated FCNC. For the case of LFV, all the 
couplings of four neutralinos will be the same order of magnitude, and all the neutralino 
contributions to LFV should be included at the same time in principle. However, we assume 
that the photino is the LSP and other neutralinos are heavy and can be ignored in order to 
simplify our analysis. It would be straighforward, although tedious, to include 4 neutralinos 
altogether and make more complete our analysis. 

There are a few differences between Ref. and our work. First of all, we can restrict 
the allowed regions of the FC mass insertion by considering different processes. Different 
processes provide independent constraints from each other, and we need not make an as- 
sumption that there is no fortuitous cancellation between 5ll and 5lr, and so on. In the 
limit the light photino dominates the LFV, we can even predict the upper bound on some 
LFV decays in a completely model independent fashion. It is straightforward to relax this 
assumption and include all the four neutralino contributions to LFV, if necessary. We con- 
sider all the LFV processes that are studied experimentally at present. We consider LFV 
decays of Z° gauge boson, and processes involving two leptons and two quarks, such as 
fi~ + Ti — > e~ + Ti, and r — > yu(e) + (a neutral meson) as well as processes involving four 
leptons and the LFV radiative decays. 

Secondly, the authors of Ref. || derived constraints on the flavor conserving mass 
insertion 5 l u from the requirement that the SUSY one-loop contribution to the lepton 
mass (one loop diagram with an insertion of 8^) is smaller than the actual lepton mass 
(Amf VSY < tn^). However, we regard this condition as an improper one, since the particle 
mass cannot be predicted by SM or SUSY models. On the contrary, it turns out that the 
anomalous magnetic moment of a lepton (ai = (g — 2)/2) can provide more meaningful and 
stronger bounds on 8 l u . 

This paper is organized as follows. In Sec. II, we construct the effective lagrangian for 
ALi = 1 and 2. The results form the basis for the calculations of transition rates for various 
LFV processes in the Sec. III. Constraints on the flavor conserving mass insertions from 
the anomalous magnetic moment are derived in Sec. IV, and the results are summarized in 
Sec. V. 



Let us first derive the effective Lagrangian for ALj = 1. A complete basis for ALi = 1 
effective Lagrangian is 



II. EFFECTIVE LAGRANGIAN FOR AL 7 = 1 AND 2 



A. C cS (4l) for ALi = 1 



AL i= l 
eff 



(4/)= E [QO. + C'M 



(6) 



i=3,5,7 



where 
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°7 = mlj^PRkF^. (7) 

The operators O^s and the associated Wilson coefficients C'^s are obtained from Oj's and 
Cj's by the exchange L ^ R. Evaluating the Feynman diagrams in Figs. 1 (the box 
diagrams) and 2 (the penguin diagrams), and matching the full amplitudes with those in 
the effective theory, we get 

m,~ \ J / LL 



m- l 
2a 



2 



m 



I 



C 7 — 2 

m~ 



We have neglected the final lepton mass rrij in the above expression, and x = m~ / m~ [|. As 
noted in Ref. ||, the Z— penguin contributions to /x — > 3e etc. are suppressed compared to 
the above by a factor of m 2 /Mf , and thus were safely ignored. Note that the 5^ and 8rl 
contribute only to O7, and not to 0%^, when we keep only dimension— 6 operators in our 
effective theory. 

The functions B^s (from the box diagrams, Fig. 1), Pj's (from the penguin diagrams, 
Fig. 2) are defined in Ref. |J, and shown below for completeness : 

1 + Ax — 5x 2 + Ax ln(x) + 2a; 2 ln(x) 



B 2 (x) = x 
P x (x 



P 2 (x) 



8(1 -a;) 4 
5 — Ax — x 2 + 2 ln(x) + 4x ln(x) 
2(1 -x) 4 ' 
1 - 6x + 18x 2 - 10a; 3 - 3x 4 + 12a; 3 ln(x) 
18(x- l) 5 : 
7 - 18a; + 9x 2 + 2a; 3 + 3 ln(x) - 9x 2 \n(x) 



9(x-l) 5 
Mi (a:) = AB x {x), 

, . -1 + 9a; + 9x 2 - 17x 3 + 18a; 2 ln(a;) + 6x 3 ln(x) 

Ms(x) = 127^1^ ' (9) 



B. £ c ff(2/ - 2g) for AL { = 1 

In order to study the fi~+ Ti — > e~+ Ti, and r — > //(or e) + (neutral meson), we need 
the effective Lagrangian for li + q lj + q where q denotes a specific quark flavor. From 
Feynman diagrams analogous to Figs. 1 and 2, we obtain 



1 Strictly speaking, the photino LSP implies that x < 1. However, we consider the case x > 1 as 
well, since it would give a rough estimate of neutralino-mediated LFV's in case the photino is no 
longer the LSP. 
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= ~% {^) LL PiWffPiti E W„<? + (L^R), (10) 

Z q=u,d,s 

(5)) ll (2B 1 (x) + B 2 (x)) TrfP L k E ^iAPl - P R )q - {L R) 

q=u,d,s 

(11) 

where we assume mj = trig for simplicity 0. Again the functions Pj's and B^s are originated 
from the penguin and the box diagrams, respectively. Note that £p e n^in an d ^-bo^ 9 can 
be obtained from Eqs. (l)-(3) by replacing ef — > efe g and ef — >■ efe 2 , respectively The 
penguin contribution £penguin contains the vector quark current, and thus can contribute to 
fi~+ Ti — > e~+ Ti, and r — > + V(= p°,(f>)- On the other hand, the box contribution 
^box 2<? depends only on the axial vector quark current so that it cannot contribute to the 
aformentoned processes, but it is relevant to the process r — > + P(= 7r°, 77). One also 
has to include the operator O7 describing lj — > lj + 7 to the above effective Hamiltonian 
when calculating physical amplitude for 21 — 2q processes. 

C. £ cS for ALi = 2 

In this subsection, we derive the effective Lagrangian for ALi = 2 for completeness. This 
Lagrangian is relevant to the muonium — > antimuonium conversion, although the resulting 
effect turns out to be too small. The relevant Feynman diagrams are shown in Fig. 3. The 
results are (we fix i = l,j = 2 in this subsection) 

^ =2 = E^ ALl=2 Q t , (12) 

t=i 

where the basis operators in the effective theory are defined as 

Qi = ejaPLfx e^ a P L fi, 

Q2 = n a P R /j &fp R fji, 

Q3 = eP R /i eP R fx, 
Qa = eP L fi eP L fi, 

Q5 = eP L ^ eP R p. (13) 

By matching the full theory amplitude with the effective amplitude, one can obtain the 
Wilson coefficients as follow : 



2 In general, the function P\,B\ and B2 should be generalized as functions of two variables, 
x = m~/m~ and y = m~/m~ because of the difference between the slepton and squark masses. 



r-2l-2q 
■'-box 



-4a 2 
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ci 



AL,=2 



a 



£jALi=2 



(*„) 

Ci (with R), 



AL i= 2 



r iAL i =2 



m ? V 12 JLR 



2xf 6 (x) 



C 3 (with L++R), 



mf 
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LL 



(8[ 2 ) rr [2/ 6 (x) + 4r/ fl (x)] - 4/ 6 (x)} . (14) 



Here, the functions /^(a;) and fe(x) are defined as 

6x(l + 3x) In x — x 3 — 9x 2 + 9x + 1 



/e(x) 
/8(x) = 



3(x-lf 
6(1 + 3a;) \nx + x 3 - 9x 2 - 9a; + 17 
6(ar- l) 5 ' 



(15) 



This completes the derivation of effective Lagrangians for ALj = 1 and 2 that will be used 
in the following sections. Also, for the purpose of Z — > Zj^, we present the amplitude for 
this decay. In this case, we need the full amplitude as given in Sec. Ill E. 



III. ANALYTIC EXPRESSIONS AND NUMERICAL ANALYSES 
FOR VARIOUS LFV PROCESSES 



A., li ^ lj -\- 'y 

The amplitude for Z, — > Zj + 7* can be written as 

M(l~ - Z~ + 7*) = ee a *uj(p - q) [q 2 la (A^P L + A*P R ) + im t a a(3 q^ (A%P l + A R P R )] u t (p) 

(16) 

with 



Af = A[(with L^R), 



A R 
A 2 



c 7 



a 



47r 2 27rra~ 



= ^(with L<r+R). 
The decay rate for U — > Zj + 7 is 



r(Zi -> Zj + 7) = ^ wit 



(|^| 2 + |^| 2 ). 



(17) 



(18) 
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Note that only the transition magnetic form factors A\ ,R contribute to the on-shell photon 
emission. The off-shell photon contribution (A\ ' R form factors) is relevant to the fi — > 3e 
and Ti — > e~+ Ti. Normalizing it to the decay rate for Zj — > IjViVj, one gets 



B(Z i - > Z i + 7 ) 



a: ' 



12tt 



M 3 (x) (5; 



+ 



+ {L R) 



(19) 



One can derive the limits on 5*-'s from the experimental upper bounds listed in Table I, 
assuming there is no fortuitous cancellations among various terms, as in Ref. pj, see Table 
II. Without such assumption, one would get a band in the ((SL)ll, (SI^lr) plane (see the 
solid lines in Fig. 4). 

B. yr + Ti -> e~ + Ti 

The muon conversion to an electron on the titanium target is one of the most sensitive 
probes of LFV that may arise from physics beyond the SM. In our case, the transition 
amplitude for this process can be expressed as 



M(ji- + Ti -> e + Ti) = - e - e [gV (A%P L + AfP R 



(20) 



R 



V X 2 e 9^a9, 
q=u,d 



where Aj'^'s are defined in Eq. (12). Note that there is no box contribution to this process, 
since only the quark vector current is important for the cohenrent conversion on the Ti 
nucleus. Also as alluded before, there is no Z— penguin contribution to this process to the 
order we are working. 

The transition rate is given by 



F(jr + Ti -> e~ + Ti) 



Aa A Z 2 \F(q 2 



- -ml)\ 2 ml 



a 1/1 1 2 1 aR j^li2 



\A^-A«\ Z + \A{ 



For the Titanium target, Z = 22, A = 48, JV = 26, Z cS = 17.6 and \F(q 2 ~ -m\ 
The experimental limit on the transition rate is given by 

r(/i~ + Ti -> e~ + Ti) < 6.1 x 1(T 13 T(/i capture in Ti), 



(21) 
0.54. 

(22) 



where the muon capture rate in Ti is r(/x capture in Ti) = (2.590 ± 0.012) x 10 6 / sec. Thus 
one gets the following upper bound : 



A± - A R \ 2 + \A?-A 



R 



iL\2 



< 



4.0 x 10- n GeV 



Thus, we get 
1 



2 



(S[ 2 ) LL (-P 1 ( X ) + M 3 ( X )) + ^L(S[ 2 ) LR M : 



1 ^ 



< 3.4 x 10" 8 GeV" 



(23) 



(24) 
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and similarly for the (L <-> R) case. This is another strong constraint that is independent 
of that from fi — > e^y. 

As we noted in the previous subsection, without assuming that there is no fortuitous 
cancellations among various terms, one would get a band in the in the ((5' ■) ll, (^I^lr)- One 
can not constrain (S^ll and (S^lr independently of each other without any assumptions. 
But, combing two different experiments, as one see in Fig. 4, we obtain two different bands 
from /i — > e7 and ji~ + Ti — > e~ + Ti. Only the shaded region are allowed. 



C. L 



3lj and k 



The amplitude for l~ — > l~ljlj(= 31 j) can be calculated from the effective Lagrangians, 
Eqs. (l)-(3). It can be written as the sum of the electromagnetic penguin and the box 
contributions : 



M 



penguin 



(Pi) [?V (A^Pl + A R P R ) + m^<7 Q % (A^P l + A r PrJ J ^(p) 



X ^ Wj(p2)7aUj(P3) ~ (pi <"> P2) 

+ B R e 2 u-( Pl )-f a P R Ui(p) uj(p 2 h a PRV j (p 3 ) 

+ [Sj(pi)7 Q 'P L Ui(p) U](p2haPRVj(P3) ~ (Pi <"> P2)] 
+ fife 2 [M-(pi)7 a P i? M i (p) U](p2haPLVj(p3) - (pi <-> p 2 )] • 

Here the box form factors fi's are given by 

2a 



£1 = 



P 



= Pf(with R), 



Tim 



s(5 l jt ) LL (2B l (x) + B 2 (x)) 



B R = fi 2 (with L<-+R). 



The decay rate for Z- 



z,. z,. z+ is 



i~i~i + 

L 3 l J l J 



a 



32vr 



mj. 



A{f + |Af| 2 - 2(Af Af + A%A H * + h.c.) 
m h 14 \ 



+ (|^| 2 + |^| 2 )[^ln 



16 

y 



2m ; 



+^(l^| 2 



1^ 



fl|2\ 



km 2 +\B R \ 2 ) 



+- (Af Pf* + A\ P 2 L * + Af fif* + Af fif + /i.e.) 
(A R fif* + Af P 2 L * + A%B R * + A^B R * + h.c.) 



(25) 



(26) 



(27) 



(28) 



In case j 7^ k, one has to remove the term with Pl <-> p 2 from the above amplitude and 
divide the fif and Pf terms by a factor of 2. Then, the decay rate becomes 
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m - ljUt) = ^ < [|^| 2 + \A* | 2 - 2(AfA? + A L 2 Af* + h.c.) 
+ (|A£| 2 + L4*| 2 ) ^81n^-12^ 

+ 1 -(\B^ + \B^)+ 1 -(\B^ + \B^) 
+^ {A{B? + A? Bf* + h.c.) 
+\ (^W + A^B** + h.c.) 
A (A*B? + A5B** + h.c.) 

- {A^B^* + A%B** + h.c)] . (29) 

We calculate the branching ratios for \i — * and /i — > 3e in the allowed region shown in 
Fig. 4 for x = 0.3 0.9, and 3.0 assuming <5's are real. For x = 0.3, the decay rate for // — > 3e 
is dominated by the term which is propotional to ( | | 2 + l^l 2 ), namely // — » e7 — > 3e. So, 
there is a strong correlation between the decay rates for fj, — > ey and /j, — > 3e, see Eq. (13) 
and Fig. 5 (a). The solid line in Fig. 5 (a) denotes this correlation. For larger x, this 
correlation becomes weaker and disappears. ^From Fig. 5, we observe that the branching 
ratio for \i — > ey — > 3e is smaller than the present upper bounds for x < 1 and the region of 
high e7) and low — > e7 — > 3e) is not allowed in the models under considerations. 

Similar analyses could be done for r — > 3e and r — > 3/i. In the case of r decay, there are 
no independent experiments like as fi~ + Ti — > e~ + Ti for fj, decay at present. So, one 

could not good predictions for r decays at present. (See the discussions below Eq. (29).) 



D. r — > £1=1,2+ (neutral meson) 

In this subsection, we consider the LFV in tau decays into a lighter lepton (e or /i) plus 
a light meson such as it , rj and p°. Different decays depend on different form factors so that 
each decay mode deserve its own study. Because of the limited numbers of tau leptons that 
have been accumulated upto now, the typical upper limits on the branching ratios of LFV 
tau decays are order of ~ 10 -6 . The limits on LFV tau decays may be improved in the future 
at Tau-Charm factories or B factories. Therefore, it is important to study every possible 
LFV tau decays in various LFV models beyond the SM. In this subsection, we consider tau 
decays into a lighter lepton (e or /i) plus one light meson such as 7i°,r],p or 0. 

The amplitude for r — > £1=1,2+ (neutral pseudoscalar meson = P such as 7i ,r], etc.) can 
be derived from the effective lagrangian (6) induced by the box diagrams : 

A 2 

M(r(k,s)^k(k',s')+P(p)) = ^ fe) (2B 1 (x) + B 2 (x)) 

x l a a P L T £e 2 (P(P)\Q^Q\0) + (L^ R). (30) 
q 

Using the PCAC relations, and assuming that \r/) = \{uu + dd — 2ss)/\/6), one gets 
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( p (p)\ H e ^7a7sg|0) = i-C P f w p a , 



(31) 



with C n o = 1,C V = l/VS and f w = 93 MeV. Then the decay rate for this decay is given by 



r(r^/, + P) 



18tt mf 



—r (2B 1 + B 2 yCj 



WJ 2 + I 



(32) 



The decay rate depends on (S i3 ) l LL and (^3)^, it does not depend on Slr^l. 

The amplitudes for r — > ^=1,2+ (neutral vector meson = F such as p°,4>) can be calcu- 
lated using the effective lagrangian (1), (5) induced by the penguin diagrams : 



lor 

M(r(k, s) - k(k',s) + V(p, e*)) = CV/ y m y All a a P L r + (L ^ R) 



2 



where 



Al = (4) iL {-AW + ^ M 3 (,)) + ^ x ^ (4)^,(4 



(33) 



(34) 



and similarly for A T R . This decay is a complete analogue of fi~ + Ti — > e~ + Ti at the parton 
level. So we expect that we can constrain 5 Li and without any assumptions combining 
t — > Z + 7 and r — > Z + V in the future. When writing the amplitude in the above form, we 
have used the definition of the vector meson decay constant fy : 



E(^l W a g|0) = C v f v m v e a \ 



(35) 



with C p o = l,C+ = -1/3 and f p0 = 153 MeV, / = 237 MeV. The decay rate for r -> k + V 
is 

f v ( m T — fiy)K + m r m v ~ 2rriy) a A 



r(r -> z< + y) = 



;C£{K| 2 + I^| 2 }- (36) 



We calculate the branching ratios for r — > / + P and r — > Z + V from the constraints shown 
in Table II. The results are shown in Table III and IV. Most of these decays are expected to 
occur with small branching ratios far below the current or future experimental search limits. 
As such, it establishes the necessary amount of tau leptons in order to probe the LFV from 
nondiagonal slepton mass matrix. 



E. Z° Id 



The photino-mediated LFV can generate the LFV decays of Z bosons. The amplitude 
for Z — > Ijlj decays are given by 



M = - 



/ a 
V2tt 



'8m 2 z G F y 



(Vl + CLi) 



V2 



X 



U(P')P L V(P)) ($i)LLX%(z) ~ m RL X_ 



mi 

RL ' 



+{{vi + <n) <-> (vi-ai), P L ^P R ), 



(37) 
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where z = m^/mj and 



Xrl(z) = - l-F 5 (z) - ^F 6 (z)) , 
mj \mj mj J 

X^) = ^(^G 1 (z) + ^G 2 (z)), 
m i \ m i m I J 

yH{z) = F x {z) + ^F 2 (z) - F 3 (z) - ^ 2 F 4 (z), 
m l m l 



Y% L {z) = ^G,(z), 



(38) 



with 



vi = -- + 2sm 2 6 w , at = ~, (39) 



and mi is the mass of the heavier lepton in the final state. The functions F^z) and Gi(z) 
are defined as follows: 

1 — x 



I rl rl—x 

— ~Fi(z) = dx dy 5 — 

m~ Jo Jo xm~ + (1 — x)m~ + y{x + y — l)m z 

1 r 1 r 1 - x , x(l — x)(l — x — y) 
— K F 2 {z) = dx dy - 5 — 3 — — — , 

m~ Jo Jo (xm~ + (1 — x)m~ + y[x + y — \)m z Y 

1 r 1 x 2 

— F s( z ) = / dx n \ — r~i 2' 

m~ Jo (1 — x)m^ + xm~ 



dx- X3{1 ~ X) 



1 r 1 

m~^ A Z Jo l "((l-x)m?+xm^) 5 

1 f 1 , r 1 -* , (1 - x)(l - x - y)(l - 2y) 

—rF 5 (z) = dx dy ^— '\- -= ^ y \ 9N9 i 

m~ jo io (m^ + (1 — :r)m~ + y[x + y — \)m z ) 2 

f 1 dx C~ X dy 2x(l - x)(l - x - y) 

Jo Jo (xm~ + (1 — x)m~ + y(x + y — l)m|) 2 ' 



(40) 



and 



[' [ l - x , (l-s)(l-2y) 

— jGAz) — I dx I dy = — = — — — 

m~ Jo Jo {xm~ + (1 — x)m~ + y[x + y — \)m z y 

1 fi-*^ -2x(l-x) 



rl rl—x 

dx dy 
Jo Jo 



m~ 2 ' Jo Jo y (xm 2 + (1 — x)m~ + y(x + y — l)m|) s 
. + z-l 



l+4^-52; 2 + 4zlog2; + 22; 2 log2; 
= = M 1 (z). 



(41) 
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The branching ratio of Z — > Ulj^i processes is given by 



Br(Z -> = Br(Z -> e + e~) * 1 " 

= Br(Z -> e+e~) 



r(Z -> e+e") 
a \ 2 1 



2W 2(v? + a?) 



(2 2 
MjUfmz) - f 3 (z)) 2 + m 3 ) LR \ 2 ^fGi(z) 



(42) 



fa - «0 2 \\{#ii)RR?mz) - F,{z)f + |(^W| 2 -^G»(z) 

where B(Z° -> e+e - ) = 3.366% and 

r(Z - e+e") = ^4 • 2(t£ + a?) . 

The upper limits on the LFV Z decays are || 

B cxp (Z -> e/i) < 2.5 x 1(T 6 , 
B exp (Z -> er) < 7.3 x 1(T 6 , 

5 cxp (Z -> //r) < 1.0 x 10" 5 , (43) 

The associated constraints on 5 z 's are so loose that they are useless. Using the constraints 
obtained in the previous subsections, we find that the upper limits on the branching ratios 
for the LFV Z decays are less than lO -7 ^ 8 -* for Z — > //(e) + r, and 10~ 10 for Z — > e/i, which 
are far below the present experimental results. Any observations of LFV Z decays with 
B > 10 -7 will indicate that the source of LFV should be different from the nondiagonal 
slepton mass matrix elements. 

F. Muonium — ► antimuonium conversion 

Now let us consider the muonium — > antimuonium conversion. The current experimental 
upper limit on the transition probability in the external magnetic field B ext = 0.1 T is |TTJ 



P cxp (M — ► M) < 8.2 x lO -11 (90%C.L.). (44) 

When this process is described the following effective Lagrangian, 

G G 

£tfT M = (W)v-A &V)v-A + (^)v+A &V)v-A > ( 45 ) 

it is known that the effective couplings are constrained as 

G— < 3.0 x 10~ 3 G F , or G + _ < 2.1 x 10" 3 G F , (46) 
assuming only one of them is nonvanishing. 
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Now we can derive the muonium (M = fi + e )— > antimuonium (M = fi e + ) conversion 
induced by the ALj = 2 effective Lagrangian, (7). Note that our model predicts that 



~ ^2 fe) 2 x /eOO (or xf 6 (x)) 
< 4.6 x 10-\5[ 2 ) 2 G F , 



(47) 



because xf 6 and fe(x) are far less than one. So the effect of the nondiagonal slepton mass 
matrix on the muonium conversion is totally negligible. In other words, if one observes 
the muonium — > antimuonium conversion, it would imply that the origin of the associated 
lepton flavor violation should be something different from what we consider in this work, for 



example R— parity violation (TT[ or dilepton gauge boson, etc. fT2fl 



IV. CONSTRAINTS ON FLAVOR CONSERVING MASS INSERTION FROM 
THE LEPTON ANOMALOUS MAGNETIC MOMENTS 



In this section, we consider the limits on the flavor conserving mass insertion (5 l u ^J ^ that 
are derivable from the anomalous magnetic moments of leptons. In Ref. H, this quantity 
was constrained from the condition Amf USY < m p h ys , where the Amf USY is calculated by a 
flavor conserving mass insertion and is finite : 



Amf USY 



2a 

4:71 



, Re (4) J(») 



(48) 



where I(x) = (—1 + x — xlnx)/(l — x) 2 . The physical mass m p h ys will be given by mbare + 
8m c ,t. + Amf USY , where 5m c . t . is the mass renormalization counter term in the MSSM. 
However its finite part is arbitrary, and one has to assume that there is no large cancellation 
between it and Amp USY in order to make use of it. In other words, renormalizable couplings 
cannot be calculated from the first principle without any ambiguity. Therefore, the condition 
that Amf USY < mphy S may be a plausable assumption, but it is by no means on the firm 
ground like the constraints considered in the previous sections. On the other hand, the 
anomalous magnetic moment of a lepton is calculable in the SM and any other renormalizable 
field theories without any ambiguity. So it is meaningful to require a£ 
we adopt in the following []. 

The anomalous magnetic moment of a lepton I (= a{) is defined as 



SM + a SUSY = a ex P;Which 



CLi 



9-2 



^(0), 



where F2(q 2 ) is the magnetic form factor of a lepton : 



M{k(p, s) k{p',s) + 7(g, e)) = u(p , s') 



F 1 {q 2 )^ + iF 2 {q 



2\ ° Hu 



2m,- 



(49) 



«i(p, *)*;(?). (so) 



3 There can be also potentially important contributions from chargino-sneutrino loop (in addition 



to the heavier neutralinos-slepton loops), especially when charginos and sneutrinos are light [13]. 
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In order to derive the SUSY contribution to the anomalous magnetic moment, one cannot 
use the effective Lagrangian presented in Sec. II A, since we assumed that the final lepton 
mass is negligible when we derived Eq. (1). One has to go back to the original expression 
for the li — > I j + 7 with m; = rrij. It is straightforward to show that the flavor- conserving 
mass insertion induces 



.SUSY 



a in* 
7r m~ 



P 3 (x) + M 3 (x) (4) 



LL rrij 



(51) 



where the function Pz{x) is given by 



1 — Qx + 3x 2 + 2x 3 — Qx 2 lnx 
6(1 -x) 4 : 



(52) 



and x = m 2 /m 2 as before, and Mi j3 (x)'s are defined in Sec. II. Here we have assumed that 
($ 1 u)ll = (5 1 u)rr, and (5-Jlr = {5 1 u )rl- The first term in Eq. ([51]) arises from the slepton- 
photino loop without any mass insertion. There is also another term proportional to mf 
coming from one insertion of A'J^. However, this is suppressed compared to the above by 
an additional factor of rrii/mj, and thus can be safely neglected. From this expression, one 
can get the constraint on x and (S 1 u )lr for a given value of mj. 



The current status of the anomalous magnetic moments for e and fx are as follow [14| : 



SM 

e 

exp 



1159652156.4(1.2) x 1(T 12 
1159652188.25(4.24) x 1(T 12 

i SM = 116591711(94) x KT 11 



31.85(4.4) x 10 



-12 



1 cxp 



1165923(8.5) x KT 



a £ ew = 5.9(8.6) x 10" 



(53) 



Here, the SM predictions for electron and muon anomalous magnetic moments include the 
one loop electroweak corrections and the two loop leading log terms, as well as QED correc- 
tions including the hadronic vacuum polarization and the hadronic light-by-light scattering 
14]. We ignored the tau anomalous magnetic moment here, since the experimental value 



begins to probe the lowest order QED correction at the present. The resulting constraints 
on (5-j)iij's for i = 1,2 are shown in Table V. Comparing with the constraints obtained in 
Ref . , our constraints are more reliable and even stronger for the case of muon. 

The imaginary part of the flavor conserving mass insertion is constrained by the electric 
dipole moment (EDM) of a lepton, as discussed in Ref. 0. The bound from electron EDM 
is very strong, 



Im (5 l n ) 



LR 



fa few x 10 



-7\ 



(54) 



V. CONCLUSION 

In conclusion, we considered the LFV in general SUSY models, where the slepton mass 
matrices are not diagonal in the basis where I — / — 7 is diagonal. We worked in the 
mass insertion approximations in which (^Ij) LR ' S constitute the suitable parameters that 
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characterize the strengths of LFV. There are strong constraints on some of these parameters 
from li — »■ lj + 7 and pT + Ti — > e~ + Ti. Using these constraints, we predict the upper 
limits on other LFV processes such as li — > 3lj, r — > /i (or e)+ (neutral meson), Z — > klj^i 
and the muonium — > antimuonium conversion. 

LFV processes considered in this work are sensitive probes of the slepton mass matrices 
which are related with the SUSY breaking mechanism. Any positive LFV signal would 
herald the existence of some new physics beyond the SM, and if the predictions in this 
work are violated, then one has to think of another source of LFV other than that through 
the nondiagonal slepton mass matrices. In particular, if one imposes the constraints from 
li — > lj + 7 and \T + Ti — > e~ + Ti, then the expected ranges for other LFV processes 
are well below the current limits and the level to be achieved in the near future. If some 
LFV processes are observed at rates higher than those predicted in this work, the source of 
LFV would not be likely to be photino-mediated. For example, presence of some R— parity 
violating couplings can lead to quantitatively different predictions |TT| from those made in 
this work. 

The constraints on the flavor conserving mass insertions were derived from the anomalous 
magnetic moments of leptons. These bounds are to be considered more sensible than those 
obtained from AmsusY < m exp, since the renormalizable couplings in the renormalizable 
field theories cannot be calculated from the given lagrangian. Our constraints still imply 
that the diagonal slepton masses should be almost degenerate, especially for the first two 
generations. One has to speculate why this should be the case in general supersymmetric 
models. 

Finally, let us comment on our assumption that the LSP is a pure photino, and other 
neutralinos are heavy enough so that their effects might be ignored. In order to do more 
complete analyses for given neutralino spectra (i.e. for given Mi,M 2 ,fi and tan/?), one can 
easily include the effects of all the neutralinos in principle, and do the similar analyses as 
presented in this work. This is possible, since the neutralino spectrum is independent of the 
slepton spectra. Our approach adopted in this work can be regarded as a first step to such 
complete analyses. The qualitative features of our predictions would not change very much. 
In other words, our predictions are expected to be correct within an order of magnitude. 
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FIGURES 




(a) 



(b) 



1. 1. 
j j 



1. 1. 

i i 



Y 



X 'i 



-X- 



(c) 



(d) 



FIG. 2. Penguin diagrams for ALj = 1. Here i,j,k,l are the generation indices. 
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FIG. 4. Allowed regions in the ((5[ 2 )ll, {$i2)lr) planes for x = 0.3. The region between solid 
lines are allowed by the present experiment for fj, — > e7 and the region between dotted lines are 
allowed by the present experiment for Ti — > e~ Ti. Combining two experiments, only the 
shaded region is allowed. 



21 



B(^^3e) X 10 

0.4 r 
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(b) 



0.75 ?->. o.;''; ; ■; XX ',. ,-''o >/■■■ \ Xx , oX:v:4\ Xx ; .x x>" XX:- o 

q rVt:"-- 1 ••! iili ' r i " i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

0.5 1 1.5 2 2.5 3 3.5 4 

(c) B(/2^>ey) X 10 

FIG. 5. The branching ratios for \i — > and \i — > 3e in the allowd region by the present 
experiment for \i — ► and fi~ Ti — > e~ Ti for x = 0.3 (a), x = 0.9 (b), and x = 3 (c) assuming 
<5's are real. The solid line in (a) denotes the correlation between the decay rates for fi —>■ e~f and 
H — ► 3e. 
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TABLES 



TABLE I. Upper limits of branching ratios for LFV processes considered in the present work. 
The muon conversion rate on the Ti atom is normalized on the muon capture rate on the Ti atom. 
For the muonium conversion, see Sec. Ill F. 



Mode 



Branching ratio 



Ref. 



H->ery 
fi — > 3e 
pT + Ti — > e~ + Ti 



3.8 x KT 11 

1.0 x KT 12 

6.1 x KT 13 



muonium conversion 



G MM < 3 -° X 10~ 3 Gf 

G+f M < 2.1 x 10- 3 G F 







t — > ej 


2.7 x 10 




3.0 x 10 


t — ► 3e 


2.9 x 10 


T — > 3/i 


1.9 x 10 


r — ► fie + e~ 


1.7 x 10 


t — > efi + p~ 


1.8 x 10 


T — > e7T° 


3.7 x 10 


T — ► ^7T° 


4.0 x 10 


r — > er? 


8.2 x 10 


r -> W 


9.6 x 10 


r — > ep° 


2.0 x 10 


r -> 


6.3 x 10 


r — > ecj) 


6.9 x 10 


t ^ p4> 


7.0 x 10 



III] 

\m 
\w 

[17| 
[17| 

m 
m 
m 
m 

\w 
\m 
\m 
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TABLE II. Limits on (5\j ) ll,rr,lr,rl from U — ► lj + 7 for m ; ~ = 100 GeV and for different 
values of x assuming there is no fortuitous cancellations among various terms. 



Prorpss 


(Hnnstrai'npr] ^ 

UllU Li J. Ctil J.AV.. \A. KJ 


X 


Tiiimit.s 




fj, —* e7 


\( ,S\ T T TfTi\ 

\ \ U 12)LL,RR\ 




4.0 x 10" 


-3 






0.9 


7.6 x 10" 


-3 






3 


1.8 x 10" 


-2 




\ \ Oi 1 T D R 7" 


0.3 


2.0 x 10" 


-6 






0.9 


2.3 x 10" 


-6 






3.0 


3.8 x 10" 


-6 




1 A» „ 1 T T J") 7-, 

\\V13) LL,RR\ 


n 3 


2 5 








0.9 


4.7 








3 


11 






\(S[ 3 )lr,rl\ 


0.3 


2.0 x 10" 


-2 






0.9 


2.4 x 10" 


-2 






3.0 


4.0 x 10" 


-2 


r -► ^7 


\{$2?)lL,Rr\ 


0.3 


2.4 








0.9 


4.5 








3.0 


10 






\(S 1 2s)lr,rl\ 


0.3 


1.9 x 10" 


-2 






0.9 


2.3 x 10" 


-2 






3.0 


3.8 x 10" 


-2 



TABLE III. Upper limits for the branching ratios of r — > / + P from the constraints shown in 
Table II. 



Process 


X 


Braching Ratio 


r — > e7r 


0.3 


0.61 x 10" 12 




0.9 


0.38 x 10~ 9 




3.0 


0.24 x 10~ 8 


T — > jU7T 


0.3 


0.56 x 10" 12 




0.9 


0.34 x 10~ 9 




3.0 


0.20 x 10~ 8 


r — > ery 


0.3 


0.20 x 10" 12 




0.9 


0.13 x 10~ 9 




3.0 


0.81 x 10~ 9 


r -► W 


0.3 


0.19 x 10" 12 




0.9 


0.12 x 10~ 9 




3.0 


0.67 x 10~ 9 
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TABLE IV. Upper limits for the branching ratios of r — > / + V from the constraints shown in 
Table II. 



Process 


X 


Braching Ratio 


r — > ep° 


0.3 


0.32 x 10~ 7 




0.9 


0.42 x 10~ 7 




3.0 


0.42 x 10~ 7 




0.3 


0.30 x 10~ 7 




0.9 


0.38 x 10~ 7 




3.0 


0.34 x 10~ 7 


r — > e0 


0.3 


0.29 x 10~ 8 




0.9 


0.39 x 10~ 8 




3.0 


0.43 x 10~ 8 


r -► 


0.3 


0.26 x 10~ 8 




0.9 


0.35 x 10~ 8 




3.0 


0.36 x 10~ 8 



TABLE V. Allowed ranges for the flavor conserving mass insertion (5 1 u )lr from the anomalous 
magnetic moment of a lepton for m ; ~ = 100 GeV. We show two cases i = 1(e) and % = 2(/i) only, 
since the anomalous magnetic moment of a tau lepton is poorly measured. 



X 




(S l 22)LR 


0.3 


-3.0 x 10" 2 2.3 x 10" 2 


-6.0 x 10~ 2 ~ 1.0 x 10" 2 


0.9 


-3.6 x 10~ 2 2.7 x 10~ 2 


-7.0 x 10~ 2 ~ 1.2 x 10~ 2 


3.0 


-5.9 x 10" 2 4.5 x 10~ 2 


-0.11 ~ 2.0 x 10~ 2 
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